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We investigate with the aid of numerical renormalization group techniques the thermoelectric 
properties of a molecular quantum dot described by the negative-(7 Anderson model. We show that 
the charge Kondo effect provides a mechanism for enhanced thermoelectric power via a correlation- 
induced asymmetry in the spectral function close to the Fermi level. We show that this effect results 
in a dramatic enhancement of the Kondo-induced peak in the thermopower of negative-(7 systems 
with Seebeck coefficients exceeding K over a wide range of gate voltages. 
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Introduction. — Thermoelectric devices currently use 
bulk materials, e.g., Si-Ge, PbTe or 612X63^"-. In the fu- 
ture, devices made of nanoscale objects, such as quantum 
dots or molecules, could offer alternatives, particularly 
for low-temperature applications, such as on-chip cool- 
ing of microprocessors or low-temperature refrigeration. 
Nanoscale objects have some potential advantages over 
their bulk counterparts, for example, in scalability or in 
their high degree of tunability (e.g., via a gate voltage), 
allowing them to be operated at optimal thermoelectric 
efficiency. Molecular quantum dots, in particular, could 
be interesting to study, since a large variety of such sys- 
tems could be fabricated and investigated for interesting 
thermoelectric properties^. 

The description of electrical and thermal transport 
through quantum dots is, however, a challenging theo- 
retical task. Electrons tunneling from the leads through 
the quasi-localized levels of the dot typically experience 
a large Coulomb repulsion on the dot, giving rise to the 
spin Kondo effect^. The latter profoundly affects trans- 
port, resulting, for example, in the lifting of Coulomb 
blockade at low temperatures for a wide range of gate 
voltages and an enhanced conductance close to the uni- 
tary limit, G ~ Ga = le? jh, for symmetric coupling to 
the leads^"^". Recent experimental and theoretical work 
has also addressed the effects of Kondo correlations on 
the thermoelectric properties of such quantum dotaii"— . 
However, the Kondo-induced enhancement of the ther- 
mopower at the Kondo temperature Tk was found to be 
very smalU^, suggesting that the spin Kondo effect, in its 
simplest manifestation, is ineffective for realizing efficient 
thermoelectric devices. 

In this Rapid Communication we consider a molecu- 
lar quantum dot with an attractive onsite Coulomb in- 
teraction, ?7 < 0, described by a negative-C/ Anderson 
impurity model, Eq. ([T]) below. Such a model has been 
used to explain the dielectric properties of amorphous 



semiconductorsii, to describe highly polarized heavy 
fermion statesi^ and to investigate the noise and non- 
equilibrium transport through negative-C/ moleculesi^. 
For a molecular quantum dot, several mechanisms could 
result in t/ < 0, for example, screening by electrons 
in metallic leads can reduce an initially repulsive local 
Coulomb interaction to negative valueaii, or, a vibrating 
molecule with a local electron-phonon interaction could 
result in a net attractive Coulomb interaction^^'^'^. For 
typically used metallic electrodes, such as gold, screen- 
ing is expected to ensure the locality of the attractive 
interaction in Eq. ([Ij . 

A negative-t/ quantum dot supports a charge Kondo 
effect in which the role of spin-up and spin-down states 
in the conventional spin Kondo effect are played by the 
non-magnetic empty and doubly occupied states of the 
doti^. As in the usual spin Kondo effect, this charge 
Kondo effect results in a renormalized Fermi liquid at low 
temperatures which has important consequences for elec- 
trical and thermal transport. It is also believed to be the 
origin of superconductivity in PbTe doped with Tl, where 
the valence skipper Tl acts as a negative- f/ center—!^. 
While some aspects of the electrical transport through a 
negative-C/ molecule have been investigated^-, the most 
interesting feature of such a system, elucidated below, lies 
in its remarkable low-temperature Kondo-induced ther- 
moelectric response which, to the best of our knowledge, 
has not been previously addressed. 

Model and calculations. — Specifically, we consider a 
quantum dot described by the following two-lead Ander- 
son impurity model 

+ Y.(^Ac..da + h.c.), (1) 
where, £d is the energy of the molecular level, U < 
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is the local Coulomb interaction, cr labels the spin, and 
a = L,R labels left and right electron lead states with 
kinetic energies eka ■ The couplings of the dot to the leads 
are denoted by Ta{io) = 271 pa{uj)\ta\'^ , where Pa{i^) = 
J2k ^i'^ ~ ^ka) is the density of states of lead a. 

The linear response transport properties can be cal- 
culated from the single-particle spectral function of the 
dot A„{uj) — —lin[Gda{^ + i^)]/7r, where Gdcr{^ + i5) = 
{{da\ dX)) is the Fourier transform of the retarded single- 
particle Green function of ((IJ . The thermopower is given 
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where / is the Fermi function, e is the electronic charge, 
and T{uj) = 2TrT(uj)'^^ Aa-{uj) is the transmission func- 
tion of the dot with T(uj) = ri^Mr^M At low tem- 
perature, a Sommerfeld expansion leads to 
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where ep — is the Fermi level of the leads. In the 
absence of a magnetic field A^{uj) — Ai{uj) = A{uj) is spin 
independent. A large thermopower at low temperature 
can be achieved by either tailoring the band structure of 
the leads to give a highly asymmetric r(a;) at ep with 
a large slope r'(e_F) (Ref. or tailoring correlations to 
yield a highly asymmetric A{uj) a.t cp with a large slope 
A'{eF), or both. We concentrate on the latter which is 
robust to details of the lead density of states, and assume 
a smooth r(a;) around ep, i-c, we take r(a;) = F = 0.01 
(in units of the half bandwidth of the leads). 

The frequency and temperature dependence of A{uj, T) 
is calculated by using the numerical renormalization 
group (NRG) method^"'. Results for U /T — —8 were ob- 
tained at gate voltages — |e|yg — [£^ + 11 /2) in the range 
\Vg\ < 8F (setting e = 1). In addition, for T = 0, we 
have compared results for occupation numbers rid with 
those from functional renormalization group (fRG)^ - and 
Bethe ansatz^- techniques (see Fig. [3] below). In the fol- 
lowing, Tk = V|C^|r/4e-'^l'^l/'*r (Ref.H) denotes the rel- 
evant low energy charge Kondo scale of ([!]). Due to the 
exponential dependence on U and F, Tk can vary by or- 
ders of magnitude, e.g. for positive-?/ systems from 1 to 
200 K22. For U = -8F, we have Tk = 2.64x lO'^F < F. 

Results. — Figure [T] shows the T = spectral func- 
tion for several gate voltages. At Vg = the pseudo-spin 
states Hd — and nd = 2 are degenerate, and the spec- 
tral function is symmetric, with a Kondo resonance of 
width 0{Tk) ai lu — Q and two Hubbard satellite peaks 
at w = Ed > and uj ~ Ed + U <Q. A finite gate voltage 
Vg induces a splitting /S.E — —2Vg of the pseudo-spin 
states which is analogous to a magnetic field in the con- 
ventional spin Kondo effect, i.e. the spectral function 
becomes highly asymmetric due to the polarizing effect 
of Vg , with Ud changing substantially from its "perfectly 
screened" value of nd = 1—. This asymmetry in the 



1 

_ (b) 


1 

T/T,= 




0.1 " 

... 1 




10 - 




— ■ 100 




/ V 




/ ^ 






V /T„ = 1 

1 g K 


1 




FIG. 1: (Color online) Main panel: T — spectral function 
for U/r = —8 for different gate voltages Vg/TK- Inset (a): 
A{ijj,T = 0) near cj — 0. Inset (b): Temperature dependence 
of A(a;,r) for Vg/TK = 1. 
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FIG. 2: (Color online) Thermopower S vs temperature at 
different gate voltages Vg/TK and U /V = —8. Inset: Com- 
parison with U > Q thermopower for Vg/TK = 10. 



single-particle spectral function with a large slope at e^, 
for both spin components, is the origin of the large ther- 
mopowers to be discussed below. The analogy to the spin 
Kondo effect in a magnetic field, can be made precise for 
the case of particle-hole symmetric bands which we con- 
sider: A particle-hole transformation on the down spins 
allows the negative-?/ Anderson model in the absence of 
a local magnetic field to be mapped onto the positive-t/ 
symmetric Anderson model in a finite local magnetic field 
B — 2ed + \U\ — —2Vg'^^, thereby explaining the highly 
asymmetric spectral function of ([T]) shown in Fig. [TJ The 
polarizing effect of finite Vg ^ B is strongest at T = 
and diminishes for T ^ Tk [see Fig. (Ub)]. In terms of 
the above analogy, this corresponds to the quenching of 
the magnetization M = (rid-x — nd\)/2 at high temper- 
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atures in the corresponding positive-C/ model in a field 
B. 

Figure [2] shows the main result of this Rapid Commu- 
nication: a dramatic enhancement of the Seebeck coeffi- 
cient induced by a finite gate voltage Vg > Tpc exceed- 
ing 50 fj,V/K for Vg > 2Tk- The maximum in the ther- 
mopower occurs on a temperature scale which correlates 
with Vg and is therefore highly tunable. Corresponding 
Seebeck coefficients for [/ > in the Kondo regime are 
insignificant (see the inset of Fig. [5]). The large enhance- 
ment in iS* is due to the correlation-induced asymmetry 
in the spectral function at finite Vg. At low tempera- 
tures, explicit calculations, within Fermi-liquid theory^, 
also shed light on this enhancement. In this limit, the 
thermopower may be expressed in terms of the occupancy 
Ud of the dot as 



S{T) 



cot(7ra£(/2) 
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with 7T 1, and 7 being the linear coefficient of spe- 
cific heat of the do^ (with 7 - 1/Tk for Vg = 0). A 
finite Vg ~ Tk polarizes the charge Kondo state, lead- 
ing to 71c; 2 for Vg > 0. This enhances a nomi- 
nally small {<S^ kB/\e\) thermopower by the large factor 
cot(7rn(i/2) ^ 1. Note also, that while a finite mag- 
netic field for [/ > also leads to asymmetric spectral 
functions A^{lo) and Ai{uj) around ep, the asymmetry 
in the Kondo regime is opposite for spin up and spin 
down. Consequently, it largely cancels in the combina- 
tion Aa- entering ^ and the thermopower is not en- 
hanced. [Furthermore, for rid — 1, the factor cot(7rrtd/2) 
is very small.] Finally, note that S(T) of the negative-?/ 
model in Fig. [5] does not exhibit a sign change with in- 
creasing temperature for any finite Vg , in contrast to the 
case oi U > 0^: The sign change of the latter is due to 
a change in slope of the spectral function at the Fermi 
level, induced by a collapse of the Kondo resonance with 
increasing temperature. This cannot occur in the U < 
model, since the spectral function remains polarized by 
a finite Vg at all relevant temperatures. 

The gate voltage dependence of the thermopower and 
electrical conductance is shown in Fig. ^a) and 3(b) at 
several temperatures. Except at T < Tk, a large Seebeck 
coefficient exceeding 50^V/K can always be realized by 
a suitable choice of gate voltage. By tuning the gate 
voltage to positive or negative values about the charge 
Kondo state at V^ = one can realize the p-type or -re- 
type legs of a thermoelectric device. Note the absence 
of a Kondo plateau in G(Vg) at T < in Fig. [S^b), 
which contrasts with the U > case, and the rapid drop 
on a scale Vg ^ Tk of G(Vg) due to the suppression of 
the Kondo state by the finite gate voltage acting as a 
magnetic field in the conventional Kondo eSeciM^^^^. 
NRG results for nd{T = 0) vs V^ compare very well with 
fRG calculations at U/T = —2, —4 and with exact Bethe 
Ansatz (BA) calculations at U/T = -8 [see Fig. El^c)]. 
Figure ini^d) shows that G(T) exhibits the typical Kondo 
scaling behavior at small gate voltages Vg < Tk- 
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FIG. 3: (Color online) Gate voltage dependence of ther- 
mopower S (a) and conductance G (b) at typical temperatures 
T /Tk- Inset (c): Dot occupation number vs gate voltage 
for U/V = -2, -4, -8. fRG results for U /V = -2, -4 agree 
with NRG to less than 8% relative error, while for U/T = —8 
NRG agrees with the Bethe ansatz (Ref. [2^ very well. Inset 
(d): Temperature dependence of G at selected gate voltages 
VJTk- 
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FIG. 4: ( Color online) The power factor S^G divided by Go = 
2e^/h vs temperature and for a range of gate voltages Vg/Tx 
for U/T — —8. The locations of Tk and F are indicated by 
arrows. Inset; Upper and lower bounds for ZT, as defined in 
the text, at several gate voltages. 



The thermoelectric efficiency of a nanoscale device is 
related to its dimensionless figure of merit defined by 
ZT = PT/K, where P = S^G is the power factor, and 
K ~ Kc + A'ph the thermal conductance due to elec- 
trons (e) and phonons (ph). For metallic leads, Kc will 
give the dominant contribution to K—, while for semi- 
conducting leads, ifph will also be important. Since, to 
the best of our knowledge, no calculation of iiTph in the 
presence of Kondo correlations is available, we discuss 
the efficiency of our system in terms of the power factor 
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(T), shown in Fig. |4l and give upper and lower bound 
estimates for ZT below. The power factor is largely inde- 
pendent of details of the leads, making it a useful quan- 
tity for future comparison with experiments. It is also a 
relevant quantity for on-chip cooling of a hot source in 
microelectronics^. For each Vg the power factor exhibits 
a maximum at a temperature which is related to Vg. The 
envelope of these curves has two maxima, one at T Tk 
for Vg « 2Tk and another at high temperatures T ^ 2T 
for ~ 4r Tff. In contrast, for U > 0, the power fac- 
tor is vanishingly small in the Kondo regime, with larger 
values being obtained only at the border between mixed 
valence and Kondo regimea^^. Turning to ZT, an up- 
per bound estimate is obtained by setting Kph = 0. A 
lower bound estimate is obtained by assuming that the 
molecule is transparent to phonons. In this case, each 
phonon mode contributes the maximum ballistic ther- 
mal conductance of kq — n^ksT /'ih^^ . For three phonon 
modes we have i^Tph — 3ko resulting in a lower bound es- 
timate for ZT. Both bounds (see the inset of Fig. 2]) show 
a maximum at a temperature T that correlates with V^, 
with the upper bound exceeding 1 for Vg/Tx 3> 1 and 
T/Tk ^1. In a real device, phonons will be inelasti- 
cally scattered, e.g. by vibrational modes of the molecule, 
thereby reducing ifph below its ballistic value, especially 
at higher temperatures where anharmonic effects become 
important. Hence, our lower bound for ZT is likely too 
stringent so that a suitable choice of gate voltage could 
allow interesting values of ZT ^ 0.5 — 1 to be achieved 



at T - IOOTk. 

Conclusions. — In summary, we investigated the ther- 
moelectric properties of a negative-[/ molecular quantum 
dot exhibiting the charge Kondo effect. A small gate 
voltage Vg > Tk is found to polarize the charge on the 
dot creating a single-particle spectral function which is 
highly asymmetric about the Fermi level. This yields 
a large enhancement of the Seebeck coefficient exceed- 
ing bQ^V/K on a temperature scale comparable to Vg. 
The device is highly tunable and allows large power fac- 
tors to be achieved at virtually any temperature by a 
suitable choice of gate voltage. In addition to the above- 
mentioned possible realizations of such devices, molecular 
complexes similar to those in Ref. [U but with valence 
skipping ions'^^ such as Bi, Tl or In, acting as negative-?/ 
centers, and attached to gold leads, could be promising 
systems to look into in the future. Reducing the dimen- 
sionality of the leads, e.g. by using carbon nanotubes^, 
could further enhance the power factor^'^. 
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